Energy dependence of elliptic flow over a large pseudorapidity range in Au+Au 

collisions at RHIC 



O 

o 

(N 

o 



(N 

O 

o 
o 



X 



B.B.Backi, M.D.Baker^, M.Ballintijn^, D.S.Barton^, R.R.Betts^, A.A.Bickley^ R.BindeF, A.Budzanowski^, 

W.Busza^, A.CarroU^, Z.Chai^, M.P.Decowski", E.Garci'a^, T.Gburek^, N.Georgei^^ K.Gulbrandsen^, S.Gushue^, 

C.Halliwell^, J.Hamblen^, M.Hauer^, G.A.Heintzelman^, CHenderson"*, D.J.Hofman^, R.S.Hollis^ R.Holynski^, 

B.Holzman^ A.Iordanova'^, E. Johnson^, J.L.Kane'', J.Katzy'''^, N.Khan^, W.Kucewicz*^, P.Kulinich'', C.M.Kuo^, 

W.T.Lin^ S.Manly*, D.McLeod^ A.C.Mignerey^, R.Nouicer^e, A.Olszewski^, R.Pak^, LCPark^, 

H.Pcrnegger'', C.Reed"', L.P.Remsberg^, M.Reuter^, C. Roland"', G.Roland'', L.Rosenberg'', J.Sagerer^, P.Sarin'', 

P.Sawicki^, H.Seals^, LSedykh^, W.Skulski^, C.E.Smitli^ M.A.Stankiewicz^, P.Steinberg^, G.S.F.Stephans^ 

A.Sukhanov^, J.-L.Tang^, M.B.Tonjes'^, A.Trzupek^, C.Vale'', G.J. van Nieuwenhuizen^ , S.S.Vaurynovich^, 

R.Verdicr", G.I.Veres", E.Wenger^, F.L.H.Wolfs®, B.Wosiek^, K.Wozniak^, A.H.Wuosmaai, B.Wyslouch'' 

^ Physics Division, Argonnc National Laboratory, Argonne, IL 60439-4843 

^ Chemistry and C-A Departments, Brookhaven National Laboratory, Upton, NY 11973-5000 

■^ Institute of Nuclear Physics PAN, Krakow, Poland 

^ Laboratory for Nuclear Science, Massachusetts Institute of Technology, Cambridge, MA 02139-4307 

^ Department of Physics, National Central University, Chung-Li, Taiwan 

® Department of Physics, University of Illinois at Chicago, Chicago, IL 60607-7059 

'' Department of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742 

* Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627 

(Dated: February 7, 2008) 

This paper describes the measurement of the energy dependence of elliptic flow for charged particles 
in Au-|-Au collisions using the PHOBOS detector at the Relativistic Heavy Ion Collider (RHIC). 
Data taken at collision energies of y^s^^Jv" = 19.6, 62.4, 130 and 200 GeV are shown over a wide 
range in pseudorapidity. These results, when plotted as a function of 77' — \rj\ — yteam, scale 
with approximate linearity throughout rj' , implying no sharp changes in the dynamics of particle 
production as a function of pseudorapidity or increasing beam energy. 



PACS numbers: 25.75.-q 

The characterization of collective flow of produced par- 
ticles by their azimuthal anisotropy has proven to be one 
of the more fruitful probes of the dynamics of heavy ion 
collisions at RHIC. The elliptic flow signal (^2) at midra- 
pidity is significant and consistent with expectations from 
hydrodynamic models at low pr [ll- It has been inter- 
preted as evidence for the production of a highly ther- 
malized state, and perhaps for partonic matter ^. At 
high pt , the observed shape of elliptic flow 01 3 is consis- 
tent with calculations incorporating jet quenching [3 and 
quark coalescence p. Interestingly, the fall of V2 with 
increasing pseudorapidity (77) [3 has been less amenable 
to understanding i^- 

Given the wide range of pseudorapidity coverage and 
energies available in PHOBOS data, it is interesting to 
examine the extent to which the shape of the flow dis- 
tributions change with energy in the frame of reference 
of one of the incoming nuclei. The multiplicity distri- 
bution has been examined by PHOBOS as a function of 
v' — \'n\ ~ Vbeam (which is an approximation of the rest 
frame of one of the nuclei) and found to be energy inde- 
pendent over a wide range of 77' J3 . Data showing such 
energy independence is said to be consistent with the 
concept of "limiting fragmentation" [lOJ, which, as used 
here, may extend well beyond the region of the collision 



normally thought of as the fragmentation region. This 
work examines the degree to which the elliptic flow in 
Au-(-Au collisions at RHIC exhibits limiting fragmenta- 
tion. 

The PHOBOS detector employs silicon pad detectors 
to perform tracking, vertex detection and multiplicity 
measurements. Details of the setup and the layout of 
the silicon sensors can be found elsewhere [ll|. Detec- 
tor components relevant for this analysis include the first 
six layers of both silicon spectrometer arms, the silicon 
vertex detector (VTX) , the silicon octagonal multiplicity 
detector (OCT), three annular silicon multiplicity detec- 
tors on each side of the collision point, and two sets of 
scintillating paddle counters. 

Monte Carlo simulations of the detector performance 
were based on the Hijing event generator 12] and the 
GEANT 3.21 [l3| simulation package, folding in the sig- 
nal response for scintillator counters and silicon sensors. 

The data used in this analysis were recorded in the 
year 2000 (^5^=130 GeV), 2001 (ys^=19.6 and 200 
GeV) and 2004 {,/s^=Q2A GeV) runs. Details on the 
event selection and signal processing can be found else- 
where [3, 13 • The majority of the data was taken with 
zero magnetic field to simplify the analysis. The excep- 
tion was at y/s^ =130 GeV, where field-on data was 
included to increase statistics. The selected data at each 
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FIG. 1: The points represent the eUiptic flow, averaged over centrahty (0-40%), as a function of rj shown separately for four 
beam energies. The boxes represent systematic uncertainties at 90% C.L. (Npart) gives the average number of participants for 
each data sample. 



of the four energies correspond to the 40% most central 
events. 

The analysis presented here is very similar to that used 
in previously published PHOBOS results at y/s^=130 
GeV 7] . It is based on the anisotropy of the azimuthal 
distribution of charged particles traversing the detector. 
At the points where charged tracks pass through an ac- 
tive silicon detector, energy is deposited in the form of 
ionization. The pad where energy is deposited is said to 
be a "hit" . This analysis is based on the "subevent" tech- 
nique where one studies the correlation of hits in one part 
of the detector with the event plane angle as determined 
by hits in a different part of the detector [lj| . 

The strength of the flow is given by the n"* Fourier 
coefficient of the particle azimuthal angle distribution, 

dN 



d{4> - V'fi) 



1 + > 2v„ cos[n(0 - Vfl)], 



(1) 



where ^u is the true reaction plane angle defined by the 
impact parameter and the beam axis. This analysis was 
confined to n = 2, the so-called elliptic flow. V2 and 
"02 (our best estimate of ipB.) are calculated as in refer- 
ence p^ . 

The analysis presented here differs from that in our 
previously published elliptic flow results Q by including 
only collisions within ±10 cm of the nominal vertex posi- 
tion (along the beam axis). This constraint encompasses 
the bulk of our data at all energies and satisfies central- 
ity determination restraints at 19.6 GeV. Unfortunately, 
such collisions occur in a region of the detector that is 
quite non-uniform in azimuth. Above and below and 
to each side of the nominal collision point are holes in 
the OCT subdetector amounting to half the azimuthal 
coverage. Particles passing into the region above or be- 
low the collision point traverse the VTX detector. The 
holes to each side of the nominal collision point prevent 
shadowing of the spectrometer detectors. 

The detector hit map was symmetrized in order to use 
the subevent technique for the flow analysis. For each 
of the holes above and below the collision point, this en- 
tailed filling the hole by mapping the inner VTX layer 
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FIG. 2: The evolution of the elliptic flow in 77 bins, averaged 
over centrality (0-40%), as a function of y/s^^- The lines 
represent fits to the data. 



hits onto a virtual OCT layer. Because of the limited 
coverage of the silicon spectrometer, the holes to each 
side were filled by linearly extrapolating across the hole 
the hit density adjacent to the hole region in cf) on an 
event-by-event basis. The latter procedure yielded an 
expected loss in sensitivity, reducing the raw, measured 
flow signal by roughly 10% relative to a detector with no 
hole in that region. 

Before analyzing each event for flow, the symmetrized 
detector hit map was weighted to correct for the rela- 
tive phase space differences between detector pads due 
to geometry (acceptance weight), and the dilution of the 
asymmetries in the hit map due to the occupancy in the 
pixelized detector (occupancy weight). 

The acceptance weights were calculated in each 77 an- 
nulus through the use of individual hit weights, w°, which 
are proportional to the inverse of the average number of 
hits in each pad. These weights were determined sepa- 
rately in bins of centrality and longitudinal vertex posi- 
tion. 



The occupancy was determined on an event-by-event 
basis from the number of occupied (-/Vqcc) a-nd unoccupied 
{Nunocc) pads in smaU sections of the detector. The occu- 
pancy weight in a given section, representing the average 
number of tracks per hit pad, was determined assuming 
a Poisson statistical distribution as lla 



Occ(?7, 0) = 



M 



l-e~^' 



(2) 



where /z=ln(l + Nocc/ Nunocc) is the average number of 
tracks per pad. This occupancy was used in concert with 
the acceptance weight to produce the overall weight for 
a given hit, 



Wi = w°Occ(ryj,(/)i), 



(3) 



which was used in the determination of "02 • 

Using the weighted and symmetrized hit map, the 
resolution-corrected elliptic flow was calculated with the 
standard subevent technique used for our earlier re- 
sults Q. The subevent regions used in the event plane 
calculation were 0.1 < j??! < 3.0 for all four energies. The 
event plane resolution was calculated separately for each 
centrality bin. The resolution correction ranged from 2 
to 3 on average, with the larger correction necessary at 
19.6 GeV. For the determination of V2 in the positive 
(negative) rj region of the detector, the subevent on the 
opposite side of midrapidity was used to evaluate 'ip2 ■ 

Monte Carlo simulations showed a residual suppres- 
sion of the flow signal, thought to be dominated by back- 
ground particles carrying no flow information and the loss 
of sensitivity due to the hit map symmetrization and the 
occupancy correction algorithm. As in our earlier work, 
this suppression was corrected using simulated data by 
comparing the output resolution corrected flow signal to 
the input flow signal for many samples of simulated data 
with different shapes and magnitudes of input flow. 

Numerous sources of systematic error were investi- 
gated, including effects due to the hit definition, hit merg- 
ing, subevent definition, knowledge of the beam orbit 
relative to the detector, shape of the dN/d?7 distribu- 
tion, hole filling procedure, vertexing algorithm, trans- 
verse vertex cuts, magnetic field configuration and sup- 
pression correction determination. The effect of these 
sources depended both on r] and centrality. In general, 
the systematic error arising from each source was deter- 
mined by varying that specific aspect of the analysis (or 
several aspects in concert) within reasonable limits and 
quantifying the change in the final V2 result as a function 
of J] and centrality. The individual contributions were 
added in quadrature to derive the 90% confidence level 
error shown in the results presented here. The systematic 
uncertainty was dominated by the suppression correction 
determination. 

The fully corrected elliptic flow signal is shown for all 
four energies in Figure ^ The values shown are consis- 
tent with previous measurements where there is energy 
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FIG. 3: The elliptic flow, averaged over centrality (0-40%), 
as a function of r; — ybeam (top) and r] + ybeam (bottom) for 
each of the four energies studied in this paper. The error bars 
represent the la statistical errors only. 



and acceptance overlap [3, [la, U^ ■ The error bars repre- 
sent the Icr statistical errors and the boxes give a mea- 
sure of the systematic error for each point at 90% con- 
fidence level. The statistical errors are somewhat corre- 
lated point-to-point due to shared event plane and event 
plane resolution determinations. 

Relative to the other energies, the data at 19.6 GeV 
comprise a smaller set of events with both smaller flow 
and multiplicity. This leads to the lack of statistical 
power at 19.6 GeV apparent in Figure ^ This, in turn, 
contributes to the large systematic errors because of the 
difficulty in separating statistical and systematic effects. 

The PHOBOS 200 GeV pr-mtegrated track-based re- 
sults agree very well with the data shown in Figure ^ 
in the available range of < 77 < 2.0 ^^. Also, the 
PHOBOS track-based elliptic flow results are consistent 
with the STAR 4-particle cumulant results as a function 
of pt [l7lll8l|. This agreement, along with the fact that 
the track-based technique is expected to have a different 
(and smaller) susceptibility to non-ffow correlations, im- 
plies the hit-based results shown here do not have a sig- 



nificant contribution from non-flow correlations, at least 
in the region |?;| <2. 

All four energies in Figure ^ show a non-boost- 
invariant, roughly triangular shape peaking at midrapid- 
ity. At the lower energies the flow seems to level off 
(and maybe even rise) at high \r]\. This might be due 
to pronounced directed flow in these regions at the lower 
energies or an effect due to the presence of participant 
nucleons. At higher energies, the participants are pushed 
further out in \ri\ and the directed flow is smaller 18). 

Figure 121 shows that the magnitude of the elliptic flow 
grows logarithmically with the beam energy for differing 
regions of |r;|. The lines represent fits to the data. 
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FIG. 4: Elliptic flow, averaged over centrality (0-40%), as a 
function of 77' = I?;! — ybeam for four beam energies. The error 
bars represent the la statistical errors only. 

Figure |31 shows the results from Figure ^ plotted in 
terms of 77 ± ybeam ■ For clarity, only the statistical errors 
are shown. For this and the following plot, the highest 
I77I points in the 19.6 GeV data are not included because 
they lack significance due to large systematic errors. 

For the purpose of examining boost invariance and the 
limiting fragmentation behavior of the elliptic flow, it 
would be best to use the rapidity (y) rather than 77. Un- 
fortunately, this is not possible with the PHOBOS de- 
tector over such a large acceptance range. The effect of 



working in rj rather than y is estimated by us and oth- 
ers 19] to be small (<15%) and is not thought to change 
the qualitative features of the shapes, though deviations 
near midrapidity might occur if plotted as a function of 

y ^ ybeam • 

Figure 01 shows the elliptic flow, as seen in Figure ^ 
where data from positive and negative rj are averaged 
and plotted as a function of 77' = I77I — ybeam for all four 
energies. The four curves scale throughout the region of 
77' overlap through mid- 77 for each energy. This scaling, 
along with the fact that the shape in 77' is approximately 
linear, implies the triangular shape of V2{ri) in Figure ^ 
and the linear evolution of V2 (In ^Js^^) shown in Fig- 
ure El 

The results in Figure 01 show that elliptic flow exhibits 
limiting fragmentation in the full range of 77', reminiscent 
of what was observed in the multiplicity |9| . The degree 
to which the elliptic flow is shown to be independent of 
energy everywhere in 77' is somewhat surprising given the 
success of hydrodynamics in describing the flow at the 
higher energies in the region near mid-77. 

In summary, these results illustrate the energy- 
independence of elliptic flow in ultra-relativistic heavy 
ion collisions over a large region of 77' throughout the 
energy reach of RHIC. The degree to which the energy 
independence of the results extends to midrapidity for the 
elliptic flow is intriguing. It is difficult to reconcile this 
fact with the common assumption that the particle pro- 
duction at midrapidity differs from that in the fragmen- 
tation region, particularly at the higher energies. These 
results are not obviously compatible with the underly- 
ing (Bjorken/Feynman) picture of the nuclear collisions 
at RHIC, supported by the success of hydrodynamics at 
midrapidity. They imply the longitudinal degree of free- 
dom is not to be treated trivially in our experimental and 
theoretical efforts to understand these collisions. 
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